Spin transport in proximity induced ferromagnetic graphene 
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Ferromagnetic insulators deposited on graphene can induce ferromagnetic correlations in 
graphene. We estimate that induced exchange splittings A ~ 5 meV can be achieved by e.g. using 
the magnetic insulator EuO. We study the effect of the induced spin splittings on the graphene 
transport properties. The exchange splittings in proximity induced ferromagnetic graphene can be 
determined from the transmission resonances in the linear response conductance or, independently, 
by magnetoresistance measurements in a spin-valve device. The spin polarization of the current 
near the Dirac point increases with the length of the barrier, so that long systems are required to 
determine A experimentally. 
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I. INTRODUCTION 

The two dimensional honeycomb lattice of graphene is 
a conceptual basis for describing carbon structures such 
as fuUerenes, carbon nanotubes and individual layers of 
graphite. The fabrication of free and stable mono-layers 
of graphene a few years ago transformed this concept into 
an experimental reality that has attracted a tremendous 
interest from the research community. l&S The low energy 
excitations of charge carriers in graphene are similar to 
massless relativistic Dirac (or rather Weyl) particles. The 
Hamiltonian 



H = -ihva-V + U{r), 



(1) 



where the velocity v k, 10® m/s is the analogue of the 
Dirac electrons speed of light (in the sense of limiting ve- 
locity) in graphene and cr — (cr^;, Cj,) is a two dimensional 
vector of Pauli matrices acting on 2-spinor states related 
to the two triangular sub-lattices constituting graphene's 
honeycomb lattice. The approximate Hamiltonian ([T]) is 
valid near the Dirac points K and K' in the reciprocal 
lattice. The two inequivalent Dirac points introduce a 
two-fold valley degeneracyf? 

The carrier concentrations are typically in the range 
10^^ — lO'^^ cm~^, corresponding to a Fermi wavelength 
oi \f ~ 50 — 100 nm.^'^ The mean free path has been 
reported to be of the order /mfp ~ 0.4 |j,m.^ With im- 
proved control over the fabrication process of graphene, 
we expect to see the realization of even cleaner samples 
with longer mean free paths. 

Spintronics aim to inject, manipulate, and detect spins 
in electronic devices. Electrical spin injection in normal 
metals is routinely achieved by contacting ferromagnets 
like Fe, Ni, and Co with normal metals such as Cu and 
Al and driving a current through the system. In semi- 
conductors, electrical spin injection is more challenging 
because of the resistance mismatch between the semicon- 
ductor and possible ferromagnetic metal contacts 4 Nev- 
ertheless, spin-injection into a semiconductor is feasible 
from a conventional ferromagnet when the interface re- 
sistance between the semiconductor and the ferromagnet 
is sufficiently large, as recently demonstrated by using Fe 



Schottky contacts in Ref. ,9. Spin injection detected via 
the GMR effect in nanotubes contacted to ferromagnets 
have also been reported. 

Graphene is clearly an interesting candidate for spin- 
tronics applications since the carrier concentration can 
be controlled by gate voltages. Also, it has a very weak 
spin-orbit interaction, leading to the possibility of rela- 
tively long spin flip lengths.^"'"- ^ In a recent experiment 
on spin injection in single layer graphene the spin flip 
length is found to be 1st ~ 1 M-m at room temperature 
in dirty samples.^- Cleaner samples are expected to have 
even longer spin flip lengths. 

We explore another possibility of spin dependent trans- 
port by envisioning that graphene is put in close proxim- 
ity to a magnetic insulator. Via the magnetic proximity 
effect, exchange splittings will be induced in graphene. 
Strong proximity induced exchange splittings due to fer- 
romagnetic insulators have been observed at EuO/Al 
interfaces-iiii^ii^ The effect was attributed to the non- 
vanishing overlap between the wave functions of the lo- 
calized moments in the magnetic insulator and the itin- 
erant electrons in the metal. The electronic wave func- 
tions can be described by atomic-like wave functions at 
the surface of thin Al filmsj^^ The spatial range is sim- 
ilar for the atomic wave functions in Al and graphene, 
so we expect the overlap between the localized moments 
and itinerant electrons in graphene at EuO / graphene in- 
terfaces to induce exchange interactions comparable to 
those observed for EuO/Al. Based on the results re- 
ported in Refs. miKniii we roughly estimate that ex- 
change splittings in graphene due to the ferromagnetic 
insulator EuO could be of the order of 5 mcV (see Ap- 
pendix]^ for details). 

In this paper, we show that proximity induced split- 
tings can be observed in the tunneling conductance as- 
sociated with a tunable barrier created by the ferromag- 
netic insulator gate on top of graphene. First, for highly 
doped barriers, we demonstrate that the splitting A can 
be directly observed from the transmission resonances in 
the conductance Moreover, for low doping of the 
barrier we show that the spin polarization of the tun- 
neling current, directly related to the spin splitting A, 
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increases with increasing length of the barrier. The spin 
polarization can be studied by magnetoresistance (MR) 
measurements in a spin valve device where two magnetic 
gates are deposited in series. Such MR measurements 
could also allow to independently determine the induced 
spin splitting A. 

This paper is organized as follows: We present a model 
of a magnetic gate in Sec. |TT1 Section JTTT] reminds the 
reader of the results obtained in Refs. and [13 for the 
conductance of a square barrier in graphene. Then we 
discuss how to obtain analytical expressions for the con- 
ductance both far from and close to the Dirac point. We 
extend the spinless situation to a spin dependent bar- 
rier with an exchange splitting A between the two spin 
channels in Sec. IIVI First, we discuss the possibilities for 
extracting the splitting A directly from the conductance 
of a single highly doped barrier. Second, we study the 
dependence of the current spin polarization on the bar- 
rier height and length. Section |V] discusses the MR effect 
in a double barrier spin valve device and discusses how 
it can be used to extract A. Finally, our conclusions are 
in Sec. IVll 



II. MODEL 

A possible way to construct a ferromagnetic gate is 
to deposit a magnetic insulator, such as EuO, on top 
of a graphene sample with a metallic gate above it (see 
Fig. [T|). So far, experimental efforts have focused on de- 
positing non-magnetic gates on graphenci^^i^ The pres- 
ence of a magnetic insulator will induce an exchange split- 
ting in graphene. The normal metal gate allows to control 
the Fermi level locally, i.e. to create a tunable barrier in 
graphene. In this way, both control of the charge and 
spin carrier concentrations can be achieved. 
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FIG. 1: A ferromagnetic insulator on top of graphene induces 
an exchange splitting in graphene. A metallic gate on top of 
the insulator controls the electrostatic potential. 

We assume in this paper that the normal metal gate 
induces a sharp potential barrier below it. This is a 
reasonable assumption provided the distance d between 
the gate and the graphene layer is much shorter than 
the Fermi wavelength Xp, which is relatively long in 
graphene, \p « 50— 100 nm.— Recently, a method for 
manufacturing top gates where the distance from the 
gate to the graphene layer is of the order oi Xp has 



been demonstrated.-^S. Observation of resonance effects 
due to sharp potential steps therefore seems feasible in 
graphene. 

The exchange interaction between the localized mag- 
netic moments in the ferromagnetic insulator and the 
spins of the electrons creates an additional spin depen- 
dent offset of the barrier potential, leading to the possi- 
bility of spin dependent tunneling. We estimate in Ap- 
pendix]^ that the exchange splitting due to the magnetic 
insulator EuO can be around 5 meV. Here we assume 
that the exchange interaction is not affected by the gate 
voltage of the top metallic gate. 



III. TUNNELING PROBABILITY 

For completeness, we first review the results for tun- 
neling through a square barrier in graphene, and follow 
the derivation in Refs. and H^. We will later extend 
this discussion to a spin dependent barrier. The charge 
carriers we consider are Dirac quasi-particles, described 
by the Hamiltonian ([1]). These quasi-particles originate 
from reservoirs to the left and to the right of the ballis- 
tic graphene sample. Ep is the Fermi energy measured 
with respect to the Dirac point of the undoped graphene 
layer. At zero temperature, the transport properties are 
governed by quasi-particles that approach a square bar- 
rier of height U and length L (see Fig. [5]) at energy Ep. 
We assume ballistic transport across the barrier, and also 
that the spin flip length 4f is much longer than the other 
length scales of the problem. Given the values for ^mfp 
and ^sf reported for graphener^*^ this regime should be 
realistic. 



III 



U 



L 



FIG. 2: Square barrier of length L 

The Hamiltonian ([T|) has the following plane wave solu- 
tions in regions /, //, and /// of Fig. [21 respectivelyfii^ 



i^(iir) 



1 

-ae~ 



-I3e' 



-ikxX 



Jkyy 



(2) 



(3) 
(4) 



The momentum of the incident particle makes an angle 
9 = aictanky/kx with the x axis. The angle of refraction. 
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i.e. the corresponding angle inside the barrier, is (j) = 
a,Tcta,n qy/q^. We consider only elastic scattering at the 
interfaces and define 



,2Nl/2 



= {hv)-'\EF\ 



and 



'\Ef-U\. 



(5) 
(6) 



The parameters a — sign{Ep) and (3 = sign{Ep — U) 
determine the wave function in the corresponding regions 
as either electron like (positive sign) or hole like (nega- 
tive sign). Translational invariance in the transverse (y) 
direction implies conservation of transverse momentum: 

^ 1y ^ kp sm9 — qp sincf). (7) 

It is convenient to introduce the dimensionless variable 

Ep^U 



(8) 



as a measure of the gate voltage U. ^ = 1 corresponds to 
the case of no barrier. Throughout the rest of the paper 
we will make the substitution u = sin 9, and we recall 
that by definition akp = Ep/hv and (3qp = {Ep — U)/hv. 

Matching the wave functions at the interfaces, (x = 
0) = '0(/7)(a; = 0) and -0(/7)(a; = L) = i}(iH){x = L), and 
solving for t gives the transmission probability T = |tp 
for a given incoming angle 0>2ii2^ 



T{u) 



where 



{e - u^) (1 - u^) + u2 (1 - 0^ sin^ (q^L) 



q^L = kpL^JS, 



sin^ 0. 



, (9) 



(10) 



Both t and T are invariant under the transformation 
ky —^ —ky as a consequence of the continuity condi- 
tion dll). 

In a real device, the sample has a finite width W. The 
allowed incoming angles 9 are therefore determined by 
the channel index n, due to the quantization of the trans- 
verse modes. This quantization condition is, for the in- 
finite mass boundary condition, ky fc„ = (ri + i) /W, 
where n are integers in the range to A'^max = [kpW/n — 
1/2J, and the transverse states are superpositions of 
states with positive and negative kyM^ Provided that 
the transverse momentum is conserved across the barrier 
interfaces, Eq. ^ is valid for systems of both finite and 
infinite width 

The conductance through the barrier for each spin in- 
dependent channel is given in the Landauer-Biittiker for- 
malism as 



2 A'max 

n=0 



(11) 



where = 2 is the valley degeneracy and Tn is the trans- 
mission probability, Eq. ^ , for a given transverse chan- 
nel fc„. When the number of channels N becomes large. 



i.e. kpW ^ 1, we can replace the summation over chan- 
nels with an integration over transverse momenta, such 
that the conductance becomes 



G = 00 du T{u) = Gog 
Jo 



with Qo defined as 



^0 = 



2e2 kpW 



(12) 



(13) 



The dimensionless conductance g in (|T2|) is plotted in 
Fig. [3] as a function of the dimensionless gate voltage 
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FIG. 3: Conductance g = G/Go as a function of ^ = 
[Ef — U) / Ef normalized to one spin channel when kpL = 14. 
The solid (blue) line shows the numerical result using (|12|) 
while the dashed (red) line is computed using the approxima- 
tion ini). Go = 2e^kFWI{h-K). 

From Eq. (|10p . we see that the longitudinal momenta in 
the barrier region, qx, can be either purely real [^"^ > u^) 
or purely imaginary (^^ < u^), corresponding to prop- 
agating and evanescent modes, respectively^SS, The con- 
tribution to the conductance from the evanescent modes 
becomes dominant around ^ = 0, and the scaling of the 
conductance with length at this point resembles that of 
a diffusive system i^^i^^ For |^| < 1, the conductance (fT^ 
can be split into the contributions from propagating and 
evanescent modes: 



g= I duT{u)+ / duT{u) 

^?prop fjcvan i 



(14) 



from which it is readily seen that the evanescent modes 
dominate in the region near ^ = as long as T{u) > 
for at least some u > |^|, (see Appendix [B] for details). 
For kpL-:^! and setting ^ = in Eq. ©: 



T{u) 



1 



cosh [kpLu) 



(15) 



This corresponds to the limit A'max ^ W/L in Ref. [22. 
Upon insertion of into the integral (fT^ , we find that 
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the conductance at the Dirac point is inversely propor- 
tional to the system length: 



(16) 



kpL 



This corresponds to the so-called minimal conductivity 
gsG X L/W — 5s5t)e^/(/i7r),^^ ffs = 2 being the spin de- 
generacy. 

For 1^1 < 1 and kpL ^ 1 we can approximate the 
conductance by the expression 

g~iai+ a^mi + exp {-kprn), (17) 

KpL 

with ai = 0.79 and 02 — 0.21 (see Appendix [Cl for details 
and Fig. [3] for a comparison with the exact solution). 
Equation (fT7)) reduces to (fTBjl when ^ 0. 

For 1^1 > 1, corresponding to a well or a large bar- 
rier, only propagating modes contribute, and we would 
expect to see resonances in the conductance due to quasi- 
bounc^ states in the barrier region. In the limit |^| ^ 1, 
using that < 1, the tunneling probability © becomes 



Tin) 



1 — cos^ {kpLC) ' 
resulting in the expression 

I cos(fcFiOI ^ siv? (kpLS^) arctanh(| cos(fci?i^)|) 



(18) 



|cos(fci.L^)|3 



(19) 



for the dimensionless conductance (see Appendix |B] for 
details). The period of 5 as a function of ^ is -K/kpL. 
Also g oscillates between 2/3 and 1. The transmission 
probability analogous to ([9|) for a square barrier in a non- 
chiral two dimensional system p 



T, 



non-chiral 



4(^2 _ y2)(i _ y2) + (1 _ ^2)2 sin^ (g^L) ■ 



(20) 

also gives oscillation with the same periodicity, but in 
this case the conductance oscillates between and 1. The 
fact that the conductance given by Eq. p9| oscillates be- 
tween 2/3 and 1 is due to the perfect tunneling of carriers 
near normal incidence in graphene. Another difference 
between graphene and a non-chiral system is that the 
transmission probability of the latter, (|20p is symmet- 
ric around ^ = 0, while the transmission probability for 
graphene, (O, depends also on the sign of ^ through the 
(1 — ^)-factor in the denominator. The asymmetry for 
the case of graphene can readily be seen in Fig. [3] 



IV. SPIN DEPENDENT BARRIER 

We now turn to a situation where the two spin chan- 
nels see barriers of different heights, i.e. the bottom of the 



conduction band at the barrier is shifted differently ac- 
cording to spin. Such a shift can arise through a Zeeman 
interaction due to an in-plane magnetic field or exchange 
field. 

The exchange term A splits the system into two sepa- 
rate subsystems according to spin. For an external mag- 
netic field B, the splitting is given by A « 2iibB. We 
introduce the spin dependent variables 



Ep-U 



Ep 



± 



A 



(21) 



where ± denotes spins parallel (-I-) or anti parallel (— ) to 
the exchange field (see Fig. [4|). In the following we will 
let g~^^~^ denote the spin resolved conductance for spins 
parallel (anti-parallel) to the exchange field. Assuming 
no spin flip, ^sf S> i, the total conductance gx across the 
barrier is given by the sum: 



(22) 
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FIG. 4; Ferromagnetic proximity effect splits the barrier ac- 
cording to spin such that = C/ =F A. 

Because A/_B « 5.8 x 10^^ nieV/T, a direct interac- 
tion of the spins with an external magnetic field gives 
only a very weak effect (about 1 meV per 20 T) , and one 
will have to rely on more indirect effects to observe such 
spin splittings. 

We propose depositing a ferromagnetic insulator, e.g. 
EuO, on top of the graphene sample to induce an ex- 
change splitting in graphene. A normal gate on top of the 
insulator allows to control the Fermi level in the same re- 
gion. The resulting potential profile is sketched in Fig. U) 
A rough estimate suggests that the splitting energy can 
be of order A w 5 meV at EuO /graphene interfaces (see 
Appendix rS| . 

As can be seen from Fig. [5] the effect of the splitting is 
simply to shift the conductance of each spin channel with 
respect to the other, leading to a broadening of the dip 
in the total conductance gx near the Dirac point ^ = 0. 
To be able to observe the splitting directly in the gx near 
the Dirac point, the magnitude of the splitting must be 
larger than the width of the dip of each spin resolved 
conductance, g^'^~^ . A measure w = (kpL)^^ of the 
width of the dip is discussed in Appendix [Cl leading to 
the condition 



L > 



hv 



(23) 



for observation of the splitting directly in gx at the Dirac 
point. However, the broadening of the dip due to a spin 
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FIG. 5: Spin resolved conductance through a square barrier 
for kpL = 14 and S = A/Ep = 0.05. The normalization of 
conductance is chosen as in Fig. [3] to correspond to g{l) — 1 
for each spin channel. 



splitting would be difficult to distinguish from a broad- 
ening due to other effects. 

From Fig. [5] it is apparent that the spin splitting has 
a more dramatic effect on the total conductance gx at 
large barrier doping, since due to the transmission reso- 
nances, and g~ can differ substantially at a given ^. 
The asymptotic expression for |^| 1 implies that 
gx has periodicity ir/kpL in ^ for (5 = 0, as shown at 
the bottom of Fig. [6l With increasing S, each peak of gx 
gradually splits into two spin resolved peaks. The split- 
ting measured from the conductance 2S equals 2A/Ef 
(see Fig.|5]), so in principle A can be determined directly 
from the total conductance across the barrier in this way. 




FIG. 6: Total conductance gr ~ g'^ + g~ when |^| 3> 1 for a 
range of different splittings S = A/Ep- For clarity the curves 
are shifted upwards in steps of 0.5 with increasing 5. 

On the other hand, it is also possible to study the 
splitting by examining the spin polarization across the 
barrier. 

We define a normalized spin polarization p along the 
direction of the exchange field as 



Inserting the approximate expression for the conductance 
from Eq. (jl7p and comparing to exact numerical calcula- 
tions, we find good agreement in the whole region I^CI < 1 
(see Fig. El). 



0.8 
0.6 
0.4 
0.2 



-0.2 
-0.4 
-0.6 
-0.8 



1 1 


1 1 1 


1 1 


Numerical data 
1 Appraxii|nate expression 



-1.5 



-0.5 



0.5 



1.5 



FIG. 7: The polarization p from the approximation (|17p com- 
pared to the exact numerical result obtained directly from 
Eq. HH). Both plots are for kpL = 14 and 5 = 0.05. 

Equation (jl7p implies that the polarization becomes 
more pronounced with increasing barrier length L (see 
Fig. [5]), owing to the fact that the evanescent modes are 
increasingly suppressed as L increases. 
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(24) 



FIG. 8: Polarization p as a function of ^ = {Ef — U)/Ef for 
different barrier lengths L with 5 = 0.05. 



V. MAGNETORESISTANCE 

Placing two magnetic gates a distance D apart in 
the graphene sample is a possible way to probe the po- 
larization p in Eq. p4)) . We assume either that D is 
much larger than the mean free path Imfp (but still much 
shorter than Igf), or that the experimental setup is re- 
alized as a three-terminal experiment, where the middle 
terminal completely randomizes the momenta between 
the two barriers (see Fig. [5]). 

Assuming that no spin flip processes take place in the 
sample, the conductance for each spin channel is found 
by treating the two barriers as resistors connected in se- 
ries. Arranging the magnetizations of the ferromagnetic 
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FIG. 9: Measuring tunneling magnetoresistance by placing 
two short barriers a distance D ^ Zmfp apart. 



barriers parallel or anti-parallel to each other, gives dif- 
ferent conductances g-^f-ff and g^^, corresponding to the 
two situations in Fig. [HI respectively. We study the po- 
larization using the "pessimistic" definition of the mag- 
netoresistance: MR = (g-ff^ — g^ij,) /g-^^. For the general 
case of different left {L) and right (i?) barriers, we obtain 



MR 



4^PL9LPRgR 



{gL + gnf - [pLgL - PRgnf ' 



(25) 



assuming that the resistance of the region D between the 
barriers is negligible compared to the typical resistances 
of the barriers. For clarity we have suppressed the sub- 
script T denoting total conductance of the left (right) 
barrier: = g+ + g-^^^. 

For identical barriers, MR reduces to the simple ex- 
pression: 



MR = p2 



(26) 



The combination of Eqs. (HH), and (US]) allows to ex- 
perimentally determine A from magnetoresistance mea- 
surements. The change of sign in the polarization, shown 
in Fig. [HI is directly related to the relative shift of the 
conductances corresponding to each spin channel. The 
coefficient MR is proportional to , which produces the 
double peak structure seen in Figs. [TUl andllll The con- 
dition for observing MR effects is also given by Eq. [531 
L > hv/ \A\. However, since the MR signal is only sensi- 
tive to the spin degree of freedom, we expect MR exper- 
iments to be a more direct probe of a spin induced split- 
ting. Any broadening of the dip introduced by sources 
other than A will also be less important, since the polar- 
ization p changes sign around ^ = 0. 

For a splitting of A = 5 meV, the condition in Eq. P5)) 
gives L > 110 nm (or equivalently kpL > 20). As can be 
seen in Fig. [51 the features in the polarization becomes 



sharper when increasing the length above this value. This 
also translates into a clearer signal in the magnetoresis- 
tance, which is plotted in Fig. [TUl and [11] for barriers of 
equal and unequal heights, respectively. 
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FIG. 10: Magnetoresistance for two barriers of equal height. 
The curves are shifted upward in steps of 0.25 for clarity. 
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FIG. 11; Same as in Fig. 1101 with one barrier being lower 
than the other - = 0.1). 

Finally, even if the top gate creates a smooth tun- 
able barrier, far from the the perfectly square potential 
discussed here, magnetoresistance measurements should 
still provide an experimental demonstration of proximity 
induced ferromagnetism in graphene, as the magnetic in- 
sulator still creates a sharp splitting of the spin up and 
spin down states in the region underneath the magnetic 
insulator. The exact dependence of the polarization p on 
the splitting A may be different in this case than the one 
presented here. 



VI. CONCLUSIONS 

We suggest using magnetic insulators deposited on top 
of graphene to create ferromagnetic graphene. The ex- 
change interaction between electrons in graphene and the 
localized magnetic moments in the insulator will give rise 
to a proximity induced exchange splitting A. We have es- 
timated that the graphene exchange splitting due to the 
magnetic insulator EuO in close proximity can be around 
A = 5 meV. 
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We have studied how the conductance of a square bar- 
rier in graphene is modified by the presence of a ferromag- 
netic insulator. We show that for large barriers or deep 
wells, 1^1 ^ 1, the splitting A can be determined directly 
from the total conductance across the barrier, provided 
that the barrier is sharp enough for transmission reso- 
nances to appear. For a barrier of length L > hvj |A|, 
where v is the Fermi velocity of the charge carriers in 
graphene, A should be observable in the polarization of 
the tunneling current near the Dirac point of the barrier, 
irrespectively of whether the barrier is smooth or sharp. 

Demonstration of proximity induced ferromagnetism 
in graphene should be possible through magnetoresis- 
tance measurements both for smooth and sharp barriers. 



Note added: After completion of this work we became 
aware of a related work by Y.G. Semenov et al.fSS, where a 
similar system with a magnetic gate is considered. Their 
work discusses the possibility of a spin FET which feasi- 
bility relies on variations of the spin splitting across the 
sample of the same order of magnitude as our estimate 
for the splitting itself. 
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APPENDIX A: ESTIMATION OF EXCHANGE 
SPLITTING AT EuO/GRAPHENE INTERFACES 

Experiments on depairing at EuO/Al interfaces sug- 
gest that the superconducting quasi-particles of Al expe- 
rience an exchange field due to the Eu^^ momentsji^ii^ 
This interaction is short ranged; essentially only the near- 
est layer of Eu^^ ions contributes. It has be shown 
that the exchange interaction between Eu^^ ions and 
charge carriers can be described as a Zeeman split- 
ting 



14.15.16.17.19 



A « cJ{S^ 



(Al) 



where c is the fractional density of Eu^^ ions to that of 
itinerant electrons in Al at the interface, J is the spatial 
average of the exchange integral and (S'z) is the average 
spin of Eu^^ ions at a given temperature. 

Perpendicular to the surface of thin Al films, the 
electronic wave functions can be well approximated by 
atomic-like wave functions.—. The spatial range of an 
atomic wave function is determined by the ratio Z/n^ 
where Z is the atomic number and n is the energy level. 
Since this ratio is approximately the same for the 3s and 



3p orbitals in Al {Z/n = 13/3 « 4.3) and the 2p orbitals 
in graphene [Z/n = 6/2 = 3), we expect the overlap be- 
tween the wave functions of localized moments and itin- 
erant electrons at EuO/graphene interfaces to be compa- 
rable to those for EuO/Al. Accordingly, we assume that 
the exchange interaction between Eu^"*" ions and itinerant 
electrons to be the same at EuO/Al and EuO/graphene 
interfaces. Ref. 17 reports the value J = 15 meV for 
Eu/Al interfaces, which also agrees with the exchange 
energy /lox = 0.1 meV estimated in Ref. [l^. 

Using a nearest neighbor distance in graphene of 
1.42 we obtain for the areal density of itinerant elec- 
trons nc ~ 40 nm^^. Similarly, the areal density of 



4 nm 



Eu -ions at the surface of EuO is 'n-^^2Jr 
Together this gives c = nEu2+ /nc Ki 10~^. 

The temperature dependence of the average spin of 
Eu'^^ ions in EuO is calculated in Ref. |3ll, showing that 
3.5 > {S^) > 3 for < T < 30 K. 

Collecting all together we arrive at the estimate 



A « 5 meV 



(A2) 



for the exchange splitting in graphene due to EuO. We 
stress that this is a very rough estimate which needs to 
be tested experimentally. 



APPENDIX B: LIMITING CASES FOR THE 
CONDUCTANCE 

a. Large potential 

For large barriers or deep wells, \Ef — C/| 3> 
^ 1. The transmission probability © then becomes 



T{u) 



(Bl) 



1 — cos^ [kpL^) 
The conductance in this case is 

Jq 1 — cos^ [kpLS) 

I cos {kpLS/) \ — sin^ ikpLS) arctanh (| cos (^^^01) 



I cos (/cfL^P 



(B2) 



which oscillates between the values 2/3 and 1 with period 
■n /kpL as a function of ^. 



b. Evanescent modes 

When \Ef — U\ <^ \Ef\: the evanescent modes dom- 
inate the transport so we neglect the contribution from 
the propagating modes. Using that |^| < u for evanescent 
modes, and that |^| ^ 1, we find 



Tiu) 



cosh (kpLu) 



(B3) 
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which is vaHd for kpL ^ 1. The conductance then be- 
comes 



1 



1 



cosh [kpLu) kpL 



(B4) 



APPENDIX C: DIP OF THE TUNNELING 
CONDUCTANCE AROUND THE DIRAC POINT 

We study how the width of the dip in the conductance 
around the Dirac point scales with barrier length L. The 
contributions from both evanescent and from propagat- 
ing modes must be considered when |^| < 1. 

The conductance due to propagating modes can be 
written as 



where 



dv 



ffprop = fiOm 



1 



(CI) 



1 



sin^ (kpL£,VT^^ 



(C2) 

When kpL ^ 1, /(^) is well approximated by a linear 
curve fli + a2£, for all |^| < 1. The function /(^) deviates 
from linearity in an oscillatory fashion in a small region 
around |^| — 0, but /(^) is allways of order unity. For 
kpL 3> 1, the conductance due to propagating waves in 
the region |^| < 1 can therefore be approximated by 



fl'prop ~ (ai + ICI, 



(C3) 



where the value of the constants ai = 0.79 and 02 = 0.21 
depend weakly on kpL when kpL ^ 1 and are found by 
fitting Eq. (jC3jl to numerical calculations. 

We have not been able to obtain an analytical expres- 
sion for the contribution due to evanescent modes. How- 
ever, we note that the contribution from T(u) in Eq. ([9|) 
to evanescent modes can be well approximated by a de- 
caying exponential function. We have fitted our numer- 
ical calculations of .govan = T{u) to an exponen- 
tially decaying function of |^|: 



Ae 



(C4) 



The constant A is found to be l/Zcj-L by letting ^ ^ and 
comparing with Eq. ()16p . Numerical evidence suggest 
that B = CkpL, with C of order unity. 

We define the width w of the dip in the conductance 
at the Dirac point as w — 2|^c|, where j^d is the value of 
1^1 for which .gprop(Cc) = 5ovan(^c)- Taking advantage of 
the fact that (?cvan decays rapidly away from |^| = 0, we 
ignore the second order term in the expression for ^prop 
for the purpose of estimating the width w. We find that 



W = 2|ee| 



1 



kpL 



(C5) 



using that 2>V(l/ai) 
function. 



1, where W is the Lambert W- 
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